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SUMMARY 


A study of the subsonic flow past an infinitely long corrugated 
circular cylinder is presented to show the relation between tiro -dimensional 
and axisymmetrical flow. In fact, a solution is obtained which contains 
as limiting cases both the Prandtl-Glauert correction for two-dimensional 
flow and the Gothert correction for flow past slender bodies of revolution. 
Included in the paper are velocity -correction formulas for a cylinder 
with a single bump and for a corrugated cylinder in the presence of walls. 


INTRODUCTION 


The present paper Is concerned with the form of the compressibility 
corrections for subsonic flow which follow from the linear -perturbation 
theory. It is now well known that there are essential differences 
between the compressibility corrections for two-dimensional flow and the 
corresponding corrections for flow about slender bodies of revolution. 

This problem has been the subject of several papers by various authors. 

(See, for example, references 1 and 2.) 

The analysis presented herein shows that the relation between two- 
dimensional and axisymmetrical flow can be clearly demonstrated in the 
solution for the flow past an infinitely long corrugated cylinder.- 1 - In 
fact, a solution is obtained which contains as limiting cases both the 
Rrandtl-Glauer fc correction for two-dimensional flow and the Gothert 
correction for flow past Blender bodies of revolution. Although the 
results for these two limiting cases are already known, the result 
obtained in the present paper shows the nature of the transition from 
one limiting case to the other. The nature of this transition has been 
treated from a different point of view in reference where the bodies 
considered consisted of a family of ellipsoids ranging from the ellipsoid 
of revolution to the infinitely long elliptic cylinder. It is of 
interest that the present example is a natural extension of th9 two- 
dimensional wavy wall treated by Ackeret in a classical paper (reference 5). 

^Dr. C. C. Lin has pointed out to the author that flow past a corrugated 
circular cylinder has previously been considered by Th. von Barman 
(reference 3 ) in a different connection as an example of the 
calculation of wave drag for supersonic flow past bodies of revolution. 
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It should he mentioned that the results presented herein were 
obtained in June 19^8, -while the author was associated with the National 
Advisory Committee for Aeronautics at Langley Air Force Base, Va. 


AXISYMMETRICAL LINEiffi -PERTURBATION FLOW PAST A CORRUGATED CYLINDER 


Let u + U and v he components of fluid velocity in the axial and 
radial directions, respectively. Let $(x,r) he the pertubation velocity 
potential in terms of which 


u 


v 



( 1 ) 


where x and r are the axial and radial directions, respectively. 
The linearized differential equation for (f) is 


C 1 



( 2 ). 


where M m is the undisturbed stream Mach number. 

Let r = a + f (x) be the equation of the meridian profile of the body 
of revolution such that |f (x) |« a and |f'(x)l'«i. The boundary condition 
at the surface of the body of .revolution is then of the following form: 



a 

Uf '(x) 



( 3 ) 


Consider now the particular case in which 

f(x) = Z6 cos (*0 

where 6 is the thickness ratio q/l of the ripple (q is the amplitude 
and 2Z is the wave length of the ripple (fig. l)). An appropriate 
solution of differential equation (2) is 
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00 

< j > =/> sinryix : 


n=l 


fr) 


(5) 


where I Q and are modified- Bessel functions of order zero and 

A n and Bn are arbitrary constants. The following properties of 
I 0 and Kq are needed: 

I 0 ’(z) = Ii(x); Ko'(x) = -K x (x) 


x « 1< 


I D (x ) £ 1, Ko(x) Rj -loggX + 0.116 
I ± (x) X E^x) 3 J 


*> 


r 

I 0 (x)3 TjU) 


J 2 so. 


x » 1< 


**,(*) 3 ^(x) a 


(6) 


When the boundary of the cylinder is given by equation (4), only the first 
term of the series (equation (5)) is needed. 


For the body of revolution in an unlimited air stream, the asymptotic 
behavior of the functions I n requires that the coefficients A n vanish. 

From equations (5), (4), and (3), it follows that the form of the perturbation 
potential caused by the ripple is 



l/l-Ko 2 %^/l-Moo 2 \ 


(7) 


Equation ( 7 ) leads to the following expression for the axial velocity u 
at the surface of the body of revolution: 


u(x,a) 


rtTJS f) 

Vl-tto 2 4) 



(8) 


/ 
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The ratio of u(x,a) for compressible and for incompressible flow is then 

Un x 4) ^(4) 

U i h-bUc 2 %^ 1 -Moo 2 ^ ) 

From equations ( 9 ) and- (6) it follows, in particular, that, 
when it Jl-K? f >>r > 


^c 1 


( 10 a) 


and, when 


i£ «1, 

l 


% 


81 + 


log 4-M ^ 2 


1.03 + log e y 


(I0h) 


Equation (10a) is of the form pf the Prandtl -Glauert correction and equation 
(lOh) is of the form of the Gothert correction. The transition between the 
two forms is supplied by equation ( 9 ). Figure 2 showB the relation between 
the results of equations ( 9 ) and (10) for a given value of the undisturbed 
stream Mach number (f^j = 0.866) . 

Note that the validity of the foregoing formulas is governed by the 
following two restrictions. First, the use of the linearized boundary 

condition (equation ( 3 )) requires that 5 = 1 and, second, that 5 « — . 

~L 1 

When the second of these two restrictions is not satisfied, it is necessary 
to satisfy the boundary condition along the line r = a + rj cos ra '& ier 

than along the line r = a and, hence the velocity correction formulas 
depend on two length ratios y and 

VELOCITY COERECTION FORMULA FOR CYLINDER WITH BUMP 


From the foregoing results the corresponding results may be deduced 
for a body of revolution having a meridian profile given by an equation of the 
form 


r 


a + 


q(X) cos(Xx)dX 


(H) 


00 
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If, now, the substitutions 



IS = T) 


are made. 


equations 


(5) and (7) yield 



sin(?a;)dX, 


(12) 


(13) 


The corresponding solution for incompressible flow is 






E 0 (Xr) 

Epfra) 


sin(X?:)dX 


( 14 ) 


A comparison of equations (it) , (13), and (ll) shows that 


</> c M = ^2 


(15) 


which is in accordance with the general results for this case (reference 2 ). 

The velocity correction formula which follows from equation (13) is of the form 


u c ( x > a ) 

a ± (x,a) 


\t](X.) 






'\J l-M^a) 


cos(Xx)dX 


Xq(X) 


E 0 (Xa) 

Ej_(Xa) 


cos(Xx)dX 


(16) 


VELOCITY' CORRECTION FORMULA FOR RIPPLE IN THE PRESENCE OF TUNNEL WALTS 


If, again, a ripple of the form of equation (4) is taken with the 
boundary condition (equation (3)) at the surface of the body of revolution, 
there is now an additional condition at the boundary of a tunnel of 
radius b, namely, for r = b 



= 0 


(17) 
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The perturbation potential <f> which satisfies equations (3) and (17) 
is of the following form: 


I 


j. mb IiODj^Or) + E 1 (gh)l 0 (pr) _ , s 
* IiODKpOa) - KiOhJIiOa) Sin(ax) 


(18) 


where 



P = \fr-M * 2 i 


(19) 


When h » 00, equation (18) reduces to equation (7). The axial 

perturbation velocity u at the surface of the body of revolution 
follows from equation (18) in the form 


u(x,a) 


aOS IiCPbjEpOa) + E 1 (pb)l 0 (Pa) 
IqOb^Cpa) - K-LObJliOa) 


( 20 ) 


On the basis of equation (20), the values of — 
a b - a 

j, — ^ — > and M^, be calculated numerically. 


may, for given values of 
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